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In this work we study the production and decay of two CP-even Higgs bosons in the general CP- 
conserving two-Higgs doublet model at the LHC. We also study the limiting case of the decoupling 
scenario. For each Yukawa version of the model, we look for the region of the parameter space where 
a signal would be seen at the LHC. Taking into account theoretical and experimental constraints 
' on the two Higgs doublet model parameters, we show that the cross section can be two orders of 

Oh' magnitude above the corresponding Standard Model cross section. We have also studied in detail 

the decay profile of the Higgs bosons and showed that interesting signatures may emerge for some 
particular values of the parameters. In some scenarios, specific triple Higgs could be measured. 

I. INTRODUCTION 



> 

CERN's Large Hadron Collider (LHC) is due to start operation in 2009 with a center of mass energy 
of 10 TeV. Therefore, testing the spontaneous symmetry breaking (SSB) mechanism and the search for 
the scalar particle or particles associated to that process is most probably postponed for at least another 

\Q ■ year. On the other hand, the Tevatron [l| has recently excluded a 160-170 GeV Standard Model (SM) 

Higgs which decays to a pair of W gauge bosons. The mass region probed will obviously grow as data is 

0^ . being collected and a Higgs boson can still be discovered at the Tevatron. 

O ■ 

If a Higgs boson is found, the next task will be to understand the underlying mechanism that gives 
mass to all known particles. Characterizing the Higgs potential is the first step towards inferring the 
nature of the SSB mechanism. There are several alternatives to the SM that can be tested at the LHC. 
The most popular are the Minimal Supersymmetric Model (MSSM) or some of its variants like two Higgs 
Doublet Model (2HDM), little Higgs models and extra dimensions models among others. It is well known 
that the CP-conserving 2HDM has in its spectrum 2 CP-even scalars, h° and H°, one CP-odd scalar A 
and a pair of charged Higgs bosons, H^. Regarding the Higgs sector, the 2HDM is a clear candidate to 
be probed at the LHC because it is the minimal version of the SM that gives rise to charged Higgs bosons 
while being consistent with all available experimental data. The difference between all these models can 
in principle be investigated by examining the processes involving scalar triple and quartic vertices. 

In the SM, the potential can be studied in processes which have at least two Higgs in final state. The 
process gg — > h°h° was studied in the SM in j^, Q, in the SM with a sequential fourth generation in Q 
and in the MSSM in Q. The SM cross section is of the order of 20 fb for a Higgs mass of 115 GeV 
and it drops steeply with increasing Higgs mass as a consequence of the fast falloff of the gluon-gluon 
luminosity. A sequential fourth generation can enhance the double Higgs production rate to ~ 300 fb 
for a Higgs mass of 115 GeV. In the MSSM, the triple coupling can be probed in a wide range of Higgs 
masses whenever the resonant decay H — > h°h° is possible. The measurement will be easier for small 
tan /3, where the trilinear coupling is larger and the heavy Higgs mass is less than twice the light Higgs 
mass in major parts of the parameter space. More detailed studies at parton level were then performed 
for the SM in [ll] and for the MSSM in H [l^. Other processes such as vector boson fusion have been 
proposed for the 2HDM and shown to be rather powerful in reconstructing the triple Higgs self coupling 

The process pp — >• h°h° is also possible in the general 2HDM. In order to be able to distinguish between 
models, such as 2HDM, MSSM, sequential fourth generation and little Higgs models, it is mandatory to 
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perform a detailed study for the double Higgs production within the general 2HDM. In a previous study, 
we have considered double Higgs production at hadron colliders in the 2HDM fermiophobic limit (12| 
and we have shown that information on the triple Higgs couplings could be extracted. In this work, we 
extend the above study to the general 2HDM. We will consider all CP-even double Higgs production, 
that is, pp — > h°h ,h H ,H°H . We will present the regions of the parameter space of a general CP- 
conserving 2HDM where this process is large enough to be probed at the LHC taking into account the 
various experimental constraints on the 2HDM parameters as well as the tree-level unitarity and vacuum 
stability constraints. We will show that in some cases, even after imposing all constraints, cross sections 
are still large enough to allow for a determination of 2HDM triple Higgs couplings. In other cases only a 
deviation from the SM can be measured and other processes would be needed to complete the picture. 
We will give particular attention to the decoupling limit where the lightest Higgs mimics the SM Higgs 
and look for the possible non-decoupling effects. 

This paper is organized as follows. In the next section we review the 2HDM potential. We then present 
all experimental and theoretical constraints that bound the 2HDM parameter space. In Section HVl we 
present our numerical results. In Section [V] we discuss the final states in the different 2HDM scenarios. 
Finally, our findings are summarized and discussed in Section [VII 



II. THE CP-CONSERVING 2HDM 



The most general two-Higgs doublet model has 14 independent parameters. It violates CP explicitly 
and it can break electric charge. Even if explicit CP- violating interactions are not present in the potential, 
spontaneous CP breaking can still occur. One way to force the CP minimum not to exist is to impose the 
exact Z 2 discrete symmetry $1 — > $1, $2 - > —$2 [13| . This symmetry can be softly broken by dimension 
two terms we define as [m^ 2 ${<I > 2 + h.c.]. In this case, we can build either a CP-conserving theory or 
one that breaks CP spontaneously. In this work the vacuum structure is chosen so that the potential is 
CP-conserving. The theory build on a minimum that does not break CP-invariance nor electric charge 
was shown to be stable at tree level [l5j . Under these constraints, the most general renormalizable 
potential which is invariant under SU(2) <E> U(l) can be written as 

y($ 1; $ 2 ) = m?$l$i + m^|$2 + (m? 2 $|$2+h.c) + -Ai($ t 1 $i) 2 + iA 2 ($J$ 2 ) 2 

+ A3($J$i)($J$ 2 ) + A 4 ($ t 1 $ 2 )(^$i) + |a 5 [(^1* 2 ) 2 + h.c] , (1) 

where $i, i = 1,2 are complex SU(2) doublets with 4 degrees of freedom each and all m 2 , Xi and m\ 2 
are real which follows from the hermiticity of the potential. From the initial eight degrees of freedom, if 
the SU(2) symmetry is broken, we end up with two CP-even Higgs states usually denoted by h° and H°, 
one CP-odd state, A , two charged Higgs boson, and three Goldstone bosons. 

The potential in Eq. ([IJ has a total of 10 parameters if one includes the vacuum expectation values. 
In a CP-conserving minimum there are two minimization conditions that can be used to fix the tree-level 
value of the parameters m\ and m 2 . The combination v 2 — vf + v 2 is fixed as usual by the electroweak 
breaking scale through v 2 = (2^/2Gf)^ 1 ■ We are thus left with 7 independent parameters, namely 
(K)i=i,...,5, ^i2: an d tan/3 = v 2 /v\. Equivalently, we can take instead the set, m^o, m H o, m^o, m,H±, 
tan/3, a and m 2 2 , as the 7 independent parameters. The angle /3 is the rotation angle from the group 
eigenstates to the mass eigenstates in the CP-odd and charged sector. The angle a is the corresponding 
rotation angle for the CP-even sector. The parameter mi2 is a measure of how the discrete symmetry 
is broken. As stated previously, the potential with m\ 2 = has an exact Z 2 symmetry and is always 
CP-conserving. 

From the above potential, Eq.([T|), we can derive the triple Higgs couplings needed for the present study 
as a function of the 2HDM parameters m^o, to^o, m^o, m H ±, tan/3, a and m 2 2 . They are the same for 
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all different Yukawa versions of the 2HDM and are given by: 
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where A^ 0?l o is the triple Higgs SM coupling, e is the electric charge, mw is the W boson mass, sw is 
the sin of the Weinberg angle and s and c were used as shorthand notation for sin and cos. As we can 
see from Eqs. (0)-([5]), all triple Higgs couplings have a quadratic dependence of the physical masses of 
the fields present in the vertex and on the soft breaking term m 2 2 . These couplings also depend strongly 
on tan/3 and a. 



The Yukawa Lagrangian is a straightforward generalization of the SM one. 
the 2HDM is 



G; 



G: 



kUR 



h.c. 



Its most general form in 



(7) 



where the G\j and Gfj- are arbitrary constants, i and j are generation indices and k is the doublet number. 
However, tree-level flavor changing neutral currents (FCNCs) are severely constrained by experiment. The 
symmetry imposed to the potential can be naturally extended to the Yukawa Lagrangian to guarantee 
that FCNCs are not present. It suffices that fermions of a given electric charge couple to no more than 
one Higgs doublet [161 ] . This can be accomplished naturally by imposing on all fields appropriate discrete 
symmetries that forbid the unwanted FCNC couplings. There are essentially four ways of doing this [13] 
and so there are four variations of the model. We define as Type I the model where only the doublet 
4>2 couples to all fermions; Type II is the model where <p 2 couples to up-type quarks and 4>i couples 
to down-type quarks and leptons; in a Type III model </>2 couples to all quarks and <f>\ couples to all 
leptons; a Type IV model is instead built such that 4>2 couples to up-type quarks and to leptons and 
4>i couples to down-type quarks. Models III and IV have been explored in a number of papers [HI and 
more recently with a renewed interest in [l9j . We present all Yukawa coupling in appendix A. Regarding 
the Higgs couplings to the quarks, the models can be grouped in pairs like (I, IV) and (II, III) as the 
only difference between these two groups is in the couplings to the leptons. Moreover, the way we have 
chosen to build the Yukawa Lagrangian is such that up type quarks have the same couplings to the Higgs 
bosons in all four models. Down type quarks have different couplings in the two groups defined above, 
(I, IV) and (II, III), as can be seen in the appendix. In the large tan /3 limit, in 2HDM-(II, III), the Higgs 
coupling to a pair of down quarks is enhanced by a factor of l/cos/3 ~ tan/3 while for 2HDM-(I, IV) 
there is no such enhancement. The processes described in this work of neutral Higgs pair production are 
either one-loop and proceed via quark loops or tree-level in the case of bb fusion. We will show that the 
theoretical constraints force tan/3 to be small, let us say smaller than 5 (we will come back to this point 
later). Because the Yukawa couplings are proportional to the quark masses and some function of the 
angles, for such small values of tan/3 the bottom loop can not compete with the top loop and the results 
for gg — > h°h° proved to be almost Yukawa-model independent. We will show that the only difference 
comes from higher order corrections that enter the calculation via the width of the Higgs bosons. 



III. THEORETICAL AND EXPERIMENTAL BOUNDS 



The parameter space of the scalar potential of the 2HDM is reduced both by theoretical constraints 
as well as by the results of experimental searches. From the theoretical constraints which the 2HDM is 
subjected to, the most important are the ones that insure tree-level unitarity and vacuum stability of 
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the theory. We also force the potential to be perturbative by requiring that all quartic couplings of the 
scalar potential, Eq. ([T]), obey |Aj| < 87r for all i. For the vacuum stability conditions that ensure that 
the potential is bounded from below, we use those from [20j . which are given by 

Ai > , A 2 > , 
y%M + A 3 +min(0, A 4 - |A 5 |) > . (8) 

The most restrictive theoretical bound comes from the full set of unitarity constraints [111, [22[ established 
using the high energy approximation as well as the equivalence theorem and which can be written as 



|a±U&±|,|c ± |,|d±|,h 2 |,|/ ± |,| ffl . 2 | <8^ (9) 



with 



o± = ^|(Ai+A 2 )±^(A 1 -A 2 ) 2 + ^(2A 3 + A 4 ) 2 |, (10) 

b ± = i|(A 1 + A 2 )± v /(A 1 -A 2 )2+4Ai|, (11) 

c± = d ± = \ |(Ai + Aa) ± ^(A 1 -A 2 )2+4A2| , (12) 

ei = (A 3 + 2A 4 -3A 5 ) , e 2 = (A 3 -A 5 ), (13) 

/+ = (A 3 + 2A 4 + 3A 5 ) , /_ = (A 3 + A 5 ), (14) 

9i = 92 = (A 3 + A 4 ) . (15) 



The 2HDM parameters are also constrained by direct experimental searches and by precision experi- 
mental data. First, the extra contributions to the 5p parameter from the extra Higgs scalars [23[ should 
not exceed the current limits from precision measurements [24l |: \5p\ < 10 -3 . Such an extra contribution 
to dp vanishes in the limit m#± = m^o. To ensure that Sp will be within the allowed range we allow 
only a small splitting between m H ± and m^o. Second, it has been shown in Ref. [25| that data from 
B —> X s j impose a lower limit of m#± > 290 GeV in Models type II and III. This constraint no longer 
applies in the case of 2HDM-I and 2HDM-IV. Third, values of tan 8 smaller than m 1 are disallowed both 
by the constraints coming from Z — >■ bb and from B q B q mixing [25] for all Yukawa versions of the model. 
On the other hand, because tan/3 can not be too large due to the theoretical constraints, except in some 
very specific scenarios, there is no need to discuss experimental constraints for very large tan/3. Finally, 
searches for Higgs bosons at colliders directly constrains their masses. The constraints on neutral Higgs 
bosons depend on their decay modes. The scenario where the light CP-even Higgs h is fermiophobic, has 
been discussed in [12j. If the neutral Higgs bosons decay mainly into fermions, the OPAL, DELPHI and 
L3 collaborations have set a limit on the h° and A masses of the 2HDM .26h28j • OPAL concluded that 
the regions 1 < m h o < 55 GeV and 3 < m A a < 63 GeV are excluded at 95% CL independently of the 
values of a and tan 8 [26( . It should be noted however that this is true as long as we are away from the 
limit a — (3 sw 0. In this limit, the cross section for the process e + e~ — > Z — > ZhP vanishes in the 2HDM 
and the above bound no longer applies. L3 set a lower limit on the CP-even mass of the order 110.3 GeV 
p7l | if the Higgs-strahlung cross section is the SM one and the Higgs decays hadronically. The DELPHI 
collaboration also studied the decay h° — > A A in e + e~ — > h°Z and e + e~ — > h°A° production and large 
portions of the (m h o, m^o) plane were excluded depending on the 2HDM suppression factor that enters 
the e + e~ — > h°Z and e + e~ — > h°A° cross sections [28| with respect to the SM. In what follows, we will 
assume that all Higgs masses are above 100 GeV unless the ratio of couplings Szzyfl* Id^zh ^ s suppressed 
down to sin(/3 — a) ~ 0.1 when no bound to m^o applies. 
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IV. RESULTS 
A. The calculation 

The double CP-even Higgs production processes pp — > SiSj (Sij = h°,H°) have tree level contributions 
from qq — ► SiSj with s-channel Higgs exchange , t-channel quark exchange and one loop contributions 
from gg — > SiSj. The most relevant contribution for process qq — > SiSj is the one where q = b because the 
Higgs couples proportionally to the quark masses. The tree level diagrams for double Higgs production 
via b-quark fusion are shown in Fig. [TJ 



s, 




b b 
(a) (b) (c) 



FIG. 1: Tree-level contribution to bb — ¥ SiSj 




FIG. 2: Triangle and box contributions to gg SiSj, where S — h°, H° and F is a generic quark. 



The process gg — > SiSj occurs only at the one- loop level and the four types of diagrams that participate 
in the process are depicted in Fig. [2] The triangle diagrams, as shown in Fig. [2] (a), represent the generic 
contribution to the process with a fermion F which can be either an up or a down quark. The second 
type of diagrams are the box contributions which come in three different topologies, with virtual up and 
down quarks being exchanged. Comparing the tree-level with the loop diagrams it is easy to understand 
that no complementarity is expected from these two contributions - if a given set of parameters originate 
a small value of a(bb — > SiSj) it will forcibly lead to a small cr{gg — > SiSj) as well. We have explicitly 
checked that the tree-level contribution is always well below the loop contribution - usually less than 5 %. 
All other quark fusion processes are negligible when compared to bb — > SiSj. 

The tree and one-loop amplitudes were generated and calculated with the packages FeynArts [1^ and 
FormCalc 30]. The scalar integrals were evaluated with LoopTools [3lJ. We have included the total 
width of the scalar particle Si = h°,H° in the calculation of the corresponding amplitude. Even if the 
dominant contribution to the one-loop amplitude comes from top and bottom diagrams, in our numerical 
results we include all quarks in the loops. In models I and IV the largest contribution is always the one 
with top quarks in the loops. In the remaining two Yukawa models the situation is not as clear - for 
instance, in the triangle loop diagram with a virtual light Higgs, S — h°, when m f ~ tana tan/3mf, 
the contribution from the bottom quark triangle loop is larger than the one with the top quark loop. 
However, the same is not true for the triangle diagram with the virtual H°. We will show that because 
tan /3 has to be small, the difference between models can only arises mainly from the a angle dependence 
which is in fact a combined dependence of the two CP-even masses, mi2, tan/3 and sin a. 



1 Bose Symmetry forbids the couplings Zh°h°, Zh°H° and ZH°H°. 
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We will start our discussion with the very interesting scenario of the decoupling limit, where the light 
Higgs resembles the SM one. We will then move to the study of the more general case and try to identify 
the regions of the 2HDM parameter space that can be probed in a more general case. 



B. The decoupling limit 

The decoupling limit of the 2HDM is usually defined as the limiting case where all scalar masses except 
the lighter CP-even Higgs become infinite and the effective theory is just the SM with one doublet (see 
32j for a discussion). In this case, the CP-even h° is the lightest scalar particle while the other Higgs 
particles, H° , A and are heavy and mass-degenerate. Using purely algebraic arguments at the tree 
level, one can derive that the main consequences of the decoupling limit is that because cos(/3 — a) — > 0, 
the CP-even h° of the 2HDM and the SM Higgs hsM have similar tree level couplings to the gauge 
bosons and to the fermions [HI, [33[ . Obviously, the decoupling limit does not rigorously apply to the 
more realistic world where the particle masses are finite. Actually, one may consider a less rigorous 
scenario, labeled as the decoupling regime [32j . where the heavy Higgs particles have masses much larger 
than the Z boson mass and may escape detection in the planned experiments. 

Several studies have been carried out looking for non-decoupling effects in Higgs boson decays and 
Higgs self-interactions. Large loop effects in h° — > 77, h° — > 7Z and h° — > bb have been pointed out 
for the 2HDM j3f| [36[ and may give indirect information on Higgs masses and the involved triple Higgs 
couplings such as h°H + H~, h°H°H°, h°A°A° and h°h°h°. The non-decoup ling contributions to the 
triple Higgs self coupling h°h°h° has been investigated in the 2HDM in Ref. |33| . revealing large non- 
decoupling effects in the 2HDM. In this section we will show that the large non-decoupling effects in the 
h°h a h° coupling will modify the double Higgs pair production pp —> h a h° cross section and make it larger 
than its SM counterpart. We will focus on the scenario where all Higgs particles of the 2HDM, except 
for the lightest CP-even Higgs, are heavy and can escape detection at the early stages of next generation 
of colliders. 

It is easy to check that in the decoupling limit the relation /3 — a —> ir/2 holds; hence, the triple 
coupling h°h°h° given in Eq. ([2]) is reduced to its SM value gh°h°h = — 3m^ /w. Furthermore, we have 
gH°h°h a — ► at tree level and so the gg — > h°h° and bb —> h°h° amplitudes reduce exactly to the SM 
results. As shown in (33|, the one- loop leading contributions from all heavy Higgs boson loops and also 
from top quark loops to the effective h°h°h° coupling, can be written as 

\*'(2HDM)- 3ml ° fl I ^ (l I I m » (l I 

KUH^HDM) - - v I 1 + Wmlov2 \ l + ^ j + ^2^2 ^ + ^ ) 

+ m k fl + 4^V- \ (16) 

where M = m 12 / ^sin /3 cos /3, M $ andpi represent the mass of the H°, A or bosons and the momenta 
of the external Higgs lines, respectively. We note that in Eq. (|16[) m^o is the renormalized physical mass 
of the lightest CP-even Higgs boson h°. In our calculation of the cross section of pp —> h°h° in the 
decoupling limit, we ignore one-loop effects due to the h°bb coupling and replace the h°h°h° coupling 
by its effective coupling given in Eq. (THJ In Fig. [3] we present the double Higgs production process for 
the LHC in the decoupling limit of the 2HDM. In the left panel we show <J gg ^ h o h o for m h n = 115 as a 
function of M$ for several values of mi 2 . In the right panel we show how the effective coupling behaves 
with respect to the corresponding SM coupling when we vary M$ and mn- Let us first note that in 
the decoupling limit, the cross section depends only on these three parameters and very mildly on tan /3 
- the sensitivity of cross section to tan/3 is mild since tan/3 only enters the h°h°h° coupling through 
M 2 = mf 2 /(sin/3 cos /3). Moreover, tan/3 is constrained by perturbative unitarity to be rather moderate. 
The abrupt cuts in the plots are due perturbative unitarity and vacuum stability constraints. We are 
considering the case m\ 2 > where the non-decoupling effects are larger. It is interesting to note that in 
the decoupling limit the cross sections can be smaller than the SM one. This is because the box diagrams 
do not depend on myi neither on m$. For small m$, there is a destructive interference between boxes and 
vertices but as m$ grows the total contribution becomes dominated by the vertex diagrams which explain 
why it is enhanced for large to$. If mi2 is very large (300 GeV or more) the box contributions become 
negligible and non-decoupling effects can be seen already for small M$. For our choice of m^a = 115 
GeV, the SM cross section for LHC is small, 0.02 pb. As is clear in Fig. |3l the 2HDM contributions can 
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FIG. 3: pp — > h°h° in the decoupling limit. On the right we show the strength of the one-loop corrected h°h°h° 
coupling relative to SM, for several values of mi 2 . On the left plot we show the production cross section for 
m h o = 115 GeV as a function of M<s> for the same values of mi 2 as well as the corresponding SM cross section. 

enhance the cross section by about one order of magnitude and reach 0.4 pb for TO12 = 300 GeV. Contrary 
to what happens in the corresponding process 77 — > h°h° [Hj], we will not be able to test non-decoupling 
effects arising from the h°H + H~ vertex at one loop level because charged Higgs bosons do not couple to 
gluons. Only non-decoupling effects from higher order corrections can be probed. In the right panel of 
Fig. 02 one can see that the coupling h°h°h° is enhanced by a factor of more than eight in some cases. 
This enhancement in the h h h coupling is due to the non decoupling effect shown in Eq. (|16l) . Finally 
we remark that the bb contribution is much smaller than the gluon fusion one and has a much higher 
uncertainty coming from the PDFs. We have also checked that cc — > h°h° contribution is much smaller 
that the bb h°h° one. 



C. The general case for pp — > h h 



The allowed parameter space of the 2HDM is restricted after imposing experimental and theoretical 
constraints. Before presenting our results we will scan the parameter space imposing all theoretical 
constraints discussed before to search for allowed regions where the values of the cross section are above 
the SM ones. The most severe constraints are the 14 tree-level unitary conditions imposed on the scalar 
potential parameters. 

As explained before, the 2HDM under study has 7 independent parameters. In this particular study, 
the mass on the pseudo-scalar, m^o and the mass of the charged Higgs, m H ± , do not play any role in 
the value obtained for the cross section. The reason is that the pure scalar couplings h°h°h°, H°H°H Q 
and H°h°h°, Eqs. ([2]), ([3]) and (UJ), do not involve those parameters. The same is obviously true for 
the Yukawa couplings. The only possible contribution from those parameters occurs when one of the 
CP-even Higgs is allowed to decay to a pair of pseudo-scalars or/and to a pair of charged Higgs bosons. 
In this case, they will influence the total cross section via a change in the total width of the CP-even 
scalar in question. However, in order to check perturbative unitarity and vacuum stability constraints, 
rriAo and mn± must be specified. The most relevant parameters we have to consider are therefore m^o , 
m#o, sin a, tan/3 and above all rriyi- We have scanned the 2HDM parameter space to look for the regions 
where the cross sections are larger than the SM one. The results are shown in Figs. H] and [5j The detailed 
parton level studies presented in concluded that although hard there is a chance of measuring the 

SM triple Higgs couplings at the LHC. For some final states the measurement could only take place with 
the tenfold increase luminosity of the Super Large Hadron Collider. Although a more detailed analysis is 
needed we will take these studies as a benchmark and look for the parameter space where cross sections 
are at least larger than the SM double Higgs production. 




1 2 3 4 5 1 1.5 2 2.5 3 3.5 4 4.5 5 



tanP tanP 

FIG. 4: The allowed regions of the 2HDM where the <r(gg —¥ h°h°) is larger than the corresponding SM cross 
section (left) and the 2HDM triple coupling h°h°h° and h°h°h° normalized to SM couplings (right). We show mi2 
as a function of tan/3 for m h o = 115 GeV, m A o = 350 GeV, m H ± — 300 GeV, m H o = 250 GeV and sin a = —0.9. 




Sina 



FIG. 5: The allowed regions of the 2HDM where the a(gg —¥ h°h°) is larger than the corresponding SM cross 
section. We have chosen m h o = 115 GeV, m A o = 350 GeV, m H ± = 300 GeV and tan/3 = 1. In the left panel 
we show mi2 as a function of m H o with sina = —0.9. In the right panel we plot mi2 as a function of sina for 
m H o= 250GeV. 



In Fig. 2]we present a scan of m\2 as a function of tan/3 for m h a = 115 GeV, m A o — 350 GeV, 
m H ± = 300 GeV, m H o — 250 GeV, tan/3 = 1 and sina = —0.9. With this choice of the charged Higgs 
boson mass, the b sj constraint is fulfilled for all 2HDM models. The scan shows the allowed region of 
the 2HDM parameter space where the cross section is larger than the SM cross section. Here we see that 
for reasonable values of the other parameters tan /3 has to be below 5. We can also witness a correlation 
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between myx and tan /3: large values of one of the variables force low values of the other variable. In the 
right panel of Fig. |4l we illustrate the size of the triple Higgs couplings h°h°h° and H°h°h° normalized 
to SM one as a function of tan/3 and myi. For small tan/3 and large m\i the couplings can be larger 
than 8 x Afo^c 

In the left panel of Fig. [5] we again show the allowed region of the 2HDM parameter space where the 
cross section is larger than the SM cross section. It is clear that large values of mi 2 are preferred because 
the self-couplings have a strong dependence on this parameter. As expected, when m#o « 2m/j0 the cross 
section reaches a maximum but all values above 2m/j0 can give rise to large cross sections depending on 
the remaining parameters. From the right panel of Fig. [5j where myi as a function of sin a is shown, we 
conclude that large values of | sincv| maximize the cross section. It should be noted however that these 
plots would change if the values of the masses were modified. 

Let us move to the discussion of double Higgs production cr(gg — > h°h°) at the LHC in the framework 
of the general 2HDM. The cases of H°H° and H°h° production will be discussed at the end of this 
section. We already know that the cross section does not depend on m^o nor on rriH± except via the 
Higgs H a or h° width. We also know that, due to the theoretical constraints, tan/3 is constrained to be 
rather moderate. From various scans over the parameters space, we conclude that in order to have large 
enhancement with respect to SM cross section one needs to take small values of tan /3 together with large 
positive m\ 2 . We start by showing how the cross section pp —> h°h° behaves with the two CP-even Higgs 
masses. This is illustrated in Fig. [SJ where the cross section a(gg — > h°h°) at the LHC via gluon fusion 
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FIG. 6: The cross section for gg h°h° with m A o = 320 GeV, m H ± = 380 GeV, sin(/3 - a) = 0.1 and tan/3 = 2. 
In the left panel we present the cross section as a function of m h o for m H o = 200 GeV while in the right panel 
we can see the variation with m H o for m h o = 120 GeV. 

is shown for three different values of 77i 12 . We have chosen sin(/3 — a) = 0.1 which allow us to probe the 
very low mass region in m^o . However, we have checked that similar values for the cross section can be 
obtained for a very wide range of sin(/3 — a). As expected, the cross section grows with myi. In the left 
panel of Fig. [6] we show a(gg — > h°h°) as a function of m^o for m^o = 200 GeV. Before the threshold, 
m H o — 2m h o, the cross section is of the order of a few pb and does not depend much on the value of 
77ii2. This is because once the channel H° — > h°h° is open it becomes dominant with a branching ratio 
close to 100 %. On the other hand, there is another triangle diagram contribution, where a virtual h° is 
exchanged, with a much smaller coupling constant. Therefore, the dependence on TO12 is mild. In this 
region, there is a clear chance of measuring the triple coupling H Q h°h° as the cross section can easily be 
two orders of magnitude above the SM (note that we show the SM line for all values of the light Higgs 
mass to serve as a reference although masses below 115 GeV are already excluded). Below the threshold, 
in the region where H° — > h°h° is closed, all contributions are now on the same footing and depend 
much more on the set of parameters chosen. For large enough 777,12 the cross section is still above the SM 
one but as one can see in Fig. [B] it can also be well below the SM value for 77112 = in the large light 
Higgs mass region. In the right panel of Fig. [6J we show a plot for the same set of parameters but now 
as a function of m^o . This plot show us exactly what we have just discussed but from a different point 
of view - that of the heavy CP-even Higgs mass. Again we see that there is a limit for what the cross 
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section can grow by increasing myx and that limit is attained when Br(H° — > h°h°) ~ 100%. Finally let 
us comment that these cross sections do not vary much when we change the Yukawa couplings. There is 
only a small difference due to the Higgs widths which is however two small to be measured at the LHC. 
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FIG. 7: The cross section for gg — > h°h° (left) and the contributing scalar couplings (right) as a function of sin a for 
different values of m H o. The parameters are m h a — 110 GeV, m A o = 320 GeV, m H ± = 380 GeV, TO12 = 340GeV, 
and tan /3 = 1. 



In Fig. [7] (left panel) we show the cross section as a function of sin a for different values of mjjo , one 
below the H° — > h°h° threshold, one above but almost on the threshold and the other one above the 
threshold. In the right panel we present the value of the couplings H°h°h° and h°h Q h° normalized to 
SM one which prove to be the main contributions by far to the value of the cross section. This was 
explicitly checked by repeating the calculation with only the triangle diagrams where either a virtual h 
or a virtual H was exchanged. Combining the information from the two figures we conclude that the 
triangle diagram with the virtual H exchanged dominates for — 1 < sin a < —0.6. Therefore in this 
range, and for this set of parameters the coupling H h°h° could be measured because the cross section 
is much larger than the corresponding SM cross section. After crossing sin a = —0.6 both couplings 
grow and now the relative importance also depends on the value of m#o: if H — > h°h° is open then 
the diagram with the virtual H still dominates and there is a good chance of measuring it; below the 
threshold, and for this particular set of parameters, it is the diagram with the virtual h° exchange that 
dominates. We conclude by saying that depending on the parameters there is a good chance of measuring 
the coupling H°h Q h° or h°h°h° and that this occurs in a large portion of the parameter space. 



D. The H°H° and H°h° production modes 



Processes gg — > H H , H h differ from gg — > h h mainly due to the fact that H is heavier than h . 
This means that there is no possibility of resonant behavior - the intermediate state is always off-shell - 
and that there is a reduction of the phase space available. We may ask ourselves what would the difference 
be if the two CP-even masses were of similar magnitude. In this case, it is the relation between sin a, mn 
and tan j3 that determines the value of the cross section, that is, it is the way these three variables are 
combined that determines how large the cross section is. Because we know that tan/3 has to be small we 
have fixed tan /3 = 1 and varied the remaining two parameters for m^o = 120 GeV and m#o = 121 GeV. 
In Fig. [8] we present all three cross sections a(gg h°h°, H°h° and H°H°) as a function of sin a for 
m h o = 120 GeV, m H o = 121 GeV, m A = 200 GeV, m H ± = 300 GeV and tan/3 = 1. In the left panel 
we take m\i = 100 GeV while we choose m\i — 300 GeV for the right panel. We start by noting that 
we have checked that even if the value of m H o is shifted up several tenths of GeV the general trends 
do not change much. It is clear that in any case the cross section can reach several hundreds of fbarn. 
However, a large ran is needed because for mi2 = 100 GeV the cross section barely reaches the SM value. 
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We can see clearly two effects due to the couplings: the huge dependence on the value of sin a due to 
the Yukawa couplings and the dependence on the relation between 77742 and sin a caused by the scalar 
couplings. We conclude that if the CP-even scalar masses are similar and not too large, the production 
of these particles is heavily conditioned by the values of 777.12 and sin a. Finally we should point out that 
in this particular limit, the two CP-even Higgs bosons would be undistinguishable experimentally and 
therefore the number of events would be obtained by summing the three contributions. 



V. HIGGS SIGNATURES 



In most of the 2HDM parameter space, the dominant process in the production of a pair of neutral 
Higgs is the triangle diagram where gg — > H° — > h°h°. As the mass of the heavy Higgs grows, the triangle 
diagram gg — > h° — > h°h° starts to play an important role. Therefore searches based on these topologies 
in the SM and in the MSSM should yield similar results. Searches based on 46 final state were performed 
in 0] for the SM and in [|, 03 for the MSSM. The final state bbr + T~ was also explored in for the 
SM. The search based on rare decays, h°h° —>■ bbjj both for the SM and for the MSSM was performed 
in Q. In this section we will present the profile of the light Higgs boson in all four Yukawa type models. 
We will then comment on the results, showing that new rare decays can be used in the 2HDM case. 

We will discuss a scenario where the light Higgs boson's mass is well below w 200 GeV and all decay 
channels that have at least one other Higgs in the final state are closed. Scenarios where the light Higgs 
boson decays to a pair of charged Higgs, a pair of CP-odd Higgs or a combination of Higgs and gauge 
bosons are obviously allowed in a general 2HDM. In the decoupling limit, the lightest CP-even Higgs 
has the SM signature. This means that close to the limit (f3 — a) « tt/2 no major differences should be 
expected in the light Higgs boson profile with respect to the SM Higgs. The same is true for model type 
I. In Fig. [9] we show the profile of the light Higgs boson in Model I as a function of its mass in the left 
panel and as a function of sin a in the right panel. As in Model I only one Higgs couples to all fermions, 
the only dependence is in the mass of the light Higgs. The branching ratios do not depend on any other 
parameters of the model. It should be noted however that when sin a = 1 (the fermiophobic limit) the 
decay to gauge bosons become dominant. In this case the light Higgs will decay mainly to two photons or 
to two W bosons depending on the Higgs mass considered and this trend can be seen in the right panel 
of Fig. [5] The values presented in the right panel, sin(/3 — a) = 0.1 and tan/3 = 1 are irrelevant for Model 
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FIG. 9: Light Higgs branching ratios in model I as a function of its mass (left panel) and as a function of sin a (right 
panel). Although irrelevant in Model I we have chosen sin(/3 — a) — 0.1, tan/3 = 1, m H o = m A o = m H ± — 300 
GeV and m 12 = GeV. 



I for the channels presented. However, they are relevant for the other models and the left panel of Fig. |H] 
also shows the qualitative behavior of Models II, III and IV for the same set of parameters. Therefore, for 
the set of values presented this plot is the same in all four models. When tan/3 = 1 and at the same time 
f3 fx a the angular dependence almost cancels resulting in a SM-like scenario. As tan /3 grows this is no 
longer true: in Model II the coupling to up quarks is suppressed resulting again in a SM-like scenario with 
the bb and the t + t~ mode dominating; in Model III all couplings are now suppressed except the ones 
to the down quarks which means that the Higgs decays almost exclusively to b-pairs; finally in Model 
IV, and again for high tan/3 the Higgs decays almost exclusively to r + r _ . We should add that, moving 
away from /3 w a, as the mass of the light Higgs grows, the decays h° — > WW* and h° —> ZZ* become 
increasingly important although their branching ratios are always below the SM values. 
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FIG. 10: Light Higgs branching ratios in models II (left panel) and III (right panel) as a function of sin a for 
tan /3 = 1. Again, although irrelevant, we have chosen the other parameters as m H o = m^o = m H ± — 300 GeV 
and mi2 = GeV. 
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Most interesting is the dependence of the light Higgs branching ratio with sin a. In Fig. [10] we show the 
light Higgs branching ratios as a function of sin a for Models II (left) and III (right) and for tan/3 = 1. 
In Model II, when sin a is close to zero, the dominant decays become cc and gg. Therefore, in this region 
a decay to non b-jets can easily reach 100 % making the detection of Higgs pairs extremely hard due to 
the overwhelming QCD background. In Model III a more pleasant situation emerges. In the same sin a 
region, the decay to t + t~ dominates reaching 50% for sin a = 0. At the same time the decay channel 
ft, — > is well above the 0.1% line turning the final state t + t~ /i + \i~ into a very good candidate. 

We should point out that these final states occur in significant regions of the parameter space as can 
be seen in the figures. Finally, note that we don't show the respective plots for large tan j3 as we would 
recover the behavior described previously, that is, in Model II (left panel) the light Higgs would decay 
exclusively to bb and r + r _ and in Model III (right panel) Br(h° —> bb) ps 100%. 

In Fig. Ill] we show the light Higgs branching ratios as a function of sin a for Model IV and for tan /3 = 1 
(left panel) and tan/3 = 4 (right panel). For tan/3 = 1 the Higgs decays mainly to b pairs. Away from 
sin a — the decay to r pairs can reach interesting values and becomes the dominant decay mode when 
| sin a | fa 1. For tan/3 = 4 both modes almost share the total branching ratio between them. Once more 
we should point out that when sin a moves away from zero, not only the decay to r pairs become the 
dominant mode, but also the decay to fi pairs reach almost 0.4% for |sina| = 1. Hence, new decay 
channels can be explored in the context of 2HDM. 

The process pp — > H°H° gives rise to very similar signatures. The only disadvantage is that H° is 
in principle heavier than h° and phase space is lost. If the masses are of the same order and below let 
us say 130 GeV, than the main decays are again bb and t + t~. For larger values of the mass again the 
decays h° — > WW* and h° —> ZZ* become dominant. In this limit they differ only in the behavior with 
sin a and tan/3 in all Yukawa models except model I. Finally, process pp — > h°H°, has in principle two 
drawbacks: first the phase space is again reduced due to the presence of a heavier Higgs in the final 
state; second we can not reduce the background by asking the two invariant masses from each pair that 
originates from a Higgs bosons (for instance two b-jets) to have a similar magnitude. 



VI. DISCUSSION AND CONCLUSIONS 



In this work we have calculated the cross section for the production processes pp — > h°h°, H°H Q , H°h°. 
We have shown that gluon fusion is always the dominant contribution. Regarding the different Yukawa 
models, little difference was found in the scanned parameter space - although the cross sections presented 
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are for models II and III they are representative of all Yukawa models. 

We started by assessing the possibility of distinguishing the lightest CP-even Higgs boson from the 
SM Higgs in the decoupling limit of the 2HDM. We concluded that the non-decoupling effects entering 
the amplitude via higher order corrections to the vertex h°h°h° can be seen for very large values of the 
soft breaking parameter to 12 or for very large values of the common mass scale M$. Taking for instance 
to 12 = 300 GeV the cross section is more than one order of magnitude above the corresponding SM value. 

The general 2HDM has too many parameters to allow for a detailed analysis taking all of them into 
account. However, for the processes studied the only relevant parameters are m.12, sin a, tan/3, m^o and 
m h o . The remaining masses have only little influence via their contribution to the CP-even Higgs widths. 
We have scanned the parameter space to conclude that, as a general trend, the largest cross sections are 
obtained for large m\i and large | sina|. The parameter tan/3 is constrained to be O(l) except for very 
specific scenarios. Note however that we have showed that even for to 12 = cross sections can be 100 
times above the SM cross section when m#o > 2t71/ 1 o. The cross sections are larger for smaller values 
of the CP-even Higgs masses and can be two orders of magnitude above the SM cross section if the 
channel H a — > h°h° is open. We have shown that there are scenarios where the triple couplings could be 
measured. 

We have then moved to the study of the other production modes H°H° and H°h°. We have showed 
that all production modes can be large and that the main variation occurs with the combination of sin a 
and mi2 parameters. In fact, we have shown that when m#o w m^o the different processes differ only by 
virtue of values of the parameters sin a and 77112 chosen as tan /3 is bounded to be small. 

In the last section we have studied the possible final states for each production mode. We have shown 
that the main decays for a light CP-even Higgs are to b-quark pairs or to t + t~ . The relative importance 
of each of these decays depend very much on the values of sin a and tan /3. In 0] it was shown that the 
best channel to measure the SM triple Higgs coupling is 2b2j for a Higgs mass below 140 GeV. In all 
Yukawa versions of the 2HDM, the branching ratio of Higgs to two photons does not differ much from 
the SM. For most of the parameter space, the decay to bb is either equal or enhanced relative to the SM. 
Therefore, this channel should again be one of the most promising to measure 2HDM triple couplings. 
The channel 2bfj, + fj,~ was also studied and shown to be one order of magnitude below 2b 2j. With larger 
values for the production cross section and the enhancement of the h — > branching ratio in some 

of the parameter space of Models III and IV, this is a channel worth exploring. Finally, other channels 
like 46, 26t + t~, and 2£>/i + /i~ could be reevaluated for the low mass region when the 2HDM cross section 
is well above the SM one. For heavier Higgs we would have two distinguish two scenarios. One, where 
decay channels with at least one other Higgs in the final state are open; in this scenario a new study 
would have to be performed [37| . The other where the decays to VFVF^*^ and ZZ^ become important, 
which means a larger h° mass and a — (3 away from zero; in this scenario, the analysis will heavily depend 
on the values of the branching ratios of the decays h° — > WW^* 1 and h° — > ZZ^*\ 
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TABLE I: Coupling constants for the fermion-scalar interactions 



Appendix A: Fermion-scalar vertices 
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